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Predictions are provided for double differential cross sections and forward-backward ratios of 
D^ production in pPb (forward) and Pbp (backward) collisions at 5.02 TeV. The effect of nuclear 
corrections on the ratio of differential cross sections ratios is estimated to be ~ (10-30)% in the 
kinematically accessible region of LHCb, and interestingly this ratio is approximately fiat with 
respect to pt(D°) due to a compensation of shadowing and anti-shadowing effects arising from 
the input nuclear PDFs. In comparison to J/'i/)(— yL~yP) measurements which have already been 
performed with the available data, the cross section for D°(— K~'k^) production is expected to be 
two-orders of magnitude higher. 


INTRODUCTION 

Measurements of D hadron production in Pb-Pb [1] 
and Au-Au [2] collisions show evidence for the suppres¬ 
sion of the differential D cross section as compared to 
references pp collisions. This suppression can be success¬ 
fully described by in-medium energy loss effects expected 
in the presence of a Quark-Gluon Plasma EHZ]. However, 
a suppression of D production in heavy ion collisions is 
also expected in the absence of a hot nuclear medium 
due to cold nuclear matter (CNM) effects [5HT2] alone. 
These effects, such as the coherent, incoherent and inelas¬ 
tic scattering in nuclear matter, as well as a modification 
of the effective parton luminosities, arise as the colliding 
constituent nucleons of the heavy ion are not free. To 
interpret the suppression of D production in heavy ion 
collisions due to hot medium effects, it is necessary to 
first quantify the size of GNM effects with independent 
measurements. One way of disentangling these effects 
is to perform measurements in pA collisions, where only 
GNM effects are expected to be present. 

Measurements of J /'ip production in pPb collisions at 
5.02 TeV, which are subject to these CNM effects, have 
been performed by both ALICE and LHCb collabora¬ 
tions [laiii]. In both cases, differential measurements 
of the ratio of J/ip production in pPb (forward) and Pbp 
(backward) collisions in the nucleon-nucleon centre-of- 
mass (COM) frame are presented, an observable defined 
as 

dRjh 
dx 

For the case of the above mentioned J/ip measurements 
the observable x denotes transverse momentum {pr) and 
rapidity (y), and a sizeable suppression is observed by 
both collaborations. The data has also been compared 
to LO and NLO predictions which incorporate CNM ef¬ 
fects through a nuclear modification of the free proton 
and neutron parton distribution functions (PDFs). The 
NLO predictions [HI [TB] , which describe both pT and y 
distributions, incorporate the NLO EPS09 nuclear PDF 


(nPDF) modifications [T7] and provide a reasonable de¬ 
scription of the data — in particular the calculation 
based on a parton energy loss model m- Although this 
indicates that the EPS09 nPDFs describe the dominant 
CNM effect, it is clear that either additional effects such 
as parton energy loss or altered heavy-quark fragmen¬ 
tation functions m or refitted nPDFs with more low-x 
data (which may in fact themselves be parameterising 
such effects) are required to describe the observed J/ip 
data. Given that the size of CNM effects in pPb lead 
collisions are important for the interpretation of the ob¬ 
served D (and J/ip [H]) suppression in Pb-Pb collisions, 
they should be validated with measurements of other fi¬ 
nal states in a similar kinematic regime^. 

It has previously been shown that pQCD predic¬ 
tions [5DH2B] provide a satisfactory description of the for¬ 
ward D production as presented by the LHCb collabora¬ 
tion for 7 TeV pp collisions [30] . In recent work [28l ES] , it 
then been demonstrated how the inclusion of this data in 
a global QCD analysis of the proton provides substantial 
improvement in the description of the low-x gluon PDF. 
The main point being that, forward D production at low- 
Pt provides sensitivity to incoming partons at moderate 
(low) values of Bjoerken-xi^( 2 ) — 2 • 10“^(5 • 10“®), where 
Xi ( 2 ) is the fraction of momentum carried by the con¬ 
stituent parton of the forward (backward) travelling pro¬ 
ton. Therefore, measurements of D production probe a 
similar kinematic regime to J/if production, albeit with 
slightly larger average values for xi, X 2 , and Q^. 

As the ratio of D®(—A'“7r+) [BD] (D° will refer to the 
sum of and mesons) and J/if [SI] pro¬ 

duction cross sections measured within the LHCb fidu¬ 
cial region in pp collisions at 7 TeV is approximately 
100, and even larger for moderate pt values, a signifi¬ 
cant improvement in the statistical precision of differen- 


^ The ALICE collaboration has presented a measurement of the 
rate of D production in pPb collisions in the central region [2]. 
As compared to a pp reference, these results are consistent with 
unity within large uncertainties of about (15-20)%. 
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tial Rfb measurements in pPh collisions can be expected 
for D hadrons as compared to J/'0. Furthermore, the 
relative systematic uncertainty of the double differential 
measurements performed in pp collisions is slightly 
smaller than those for the corresponding J/'ijj measure¬ 
ment. One therefore naively expects that differential Rfb 
measurements of D hadrons will have improved precision 
as compared to the J/ip [Mj counter-part. 

In this letter NLO-I-LL predictions are provided for 
double differential D production within the LHCb fidu¬ 
cial region for both forward and backward pPb collisions 
at 5.02 TeV. Throughout, CNM effects are incorporated 
via EPS09 nPDF modifications [T7]. Finally, differential 
predictions for the observable Rfb are also provided. 


HEAVY-QUARK PRELIMINARIES 

The centrality integrated cross-section for heavy-quark 
pair production in pA collisions can be computed by ap¬ 
plying the standard factorisation formula as 

^pA^QQ ~ ^ ^ J Pp)Rjf^ 0 fj/N{Xj, 

da^j^QQX [s,pp,pp,as{pp),mQ) -I- ... 

. 

where the ellipses denote non-factorisable corrections 
which are neglected in this approximation. This differs 
from the corresponding cross section in pp collisions by 
the introduction of the term ® fj/B{xj,p%), which 
represents a flavour and mass number (A) dependent nu¬ 
clear modification factor to the parton distribution func¬ 
tions of the colliding nucleon, and a linear scaling of the 
cross section with A. 

Throughout this work, the EPS09 nPDF modifications 
are adopted as a baseline. These modification parame¬ 
ters, which are applied to free proton or neutron PDFs, 
are provided in a Hessian basis which allows the un¬ 
certainties from the global EPS09 analysis to be prop¬ 
agated to observables. The colliding Pb atom will also 
be approximated to be wholly constituted of protons, 
which is reasonable for the QCD production of heavy- 
quarks which is insensitive to the flavour of high-a; va¬ 
lence quarks. 

It is worth commenting that although the EPS09 
nPDFs have been extracted using CTEQ6.1M PDFs [32] . 
the application of these nPDFs to other input proton 
PDFs has been shown to have little to no impact on pPh 
to pp cross section ratio observables such as i?ppb [IS]- 
Checks with the LO nPDF modifications applied to 
NNPDF3.0 LO [33] and cteqhll [34] proton PDFs in¬ 
dicate this is also true for Rfb- As a baseline PDF 
set, the EPS09 nPDFs are applied to the proton VFNS 
NNPDF3.0LHCb NLO PDF set with as{mz) = 0.118 
from [2S]. This PDF set was obtained by including 


7 TeV LHCb and data [SO] into the NNPDF3.0 
global data set via the standard Bayesian reweighting 
technique [35l [36] . This PDF set has the advantage of 
an improved description of the behaviour of the gluon 
PDF at low-a;, and is also provided in the same scheme 
as the nPDFs. All PDFs are accessed via the LHAPDF6 
interface m- 

The heavy-quark predictions are computed at NLO 
using the massive calculation [38] implemented in 
POWHEG-BOX [35]. This calculation is performed in 
a fixed-flavour scheme (FFS), where the heavy-quark 
flavour is considered massive, and not included as a de¬ 
gree of freedom in the running of ag or in the PDFs. 
To consistently convolute this calculation with the five- 
flavour PDFs and nuclear corrections, it is necessary to 
perform a change of renormalisation scheme as explained 
in [25] . For charm production, it is necessary to add the 
following corrections factors in the region p > 


-( 0 ) ‘^TpUsipp) 
37r 

yol ^Tpag{pl) 



Log 



+ Log 



(3) 


where represents the Born contribution. The first 
term accounts for the change in the running of ag with 
Uf + 2 active flavours, and the second term accounts for 
both ag modifications and the depletion of the gluon PDF 
due to 5 —> QQ splittings above charm and bottom quark 
mass thresholds. In the kinematic region rric < p^ < rn^, 
only the correction factors for one flavour are neces¬ 
sary. Predictions provided in this way are consistent 
with those obtained in a FFS to sub-leading terms of 
0{a'^y. This calculation is then subsequently interfaced 
to Pythias gUSg using the POWHEG method [33533] . 
and the default PythiaS Monash 2013 tune [46] is used 
throughout. The hadronisation of heavy-quarks in this 
set-up is described by a non-perturbative model for which 
the modelling parameters have been tuned to LEP data. 
The predictions of D hadrons in this work is therefore 
accurate at NLO-I-LL accuracy, due to the collinear re¬ 
summation in the parton shower. 

When providing predictions for observables, the follow¬ 
ing sources of theoretical uncertainty will be considered: 
charm mass, nPDFs, proton PDFs, and missing higher- 
order corrections. The charm quark pole mass is taken 
to be rric = (1.5 ± 0.2) GeV and the corresponding un¬ 
certainty is found by taking the envelope of predictions 
found with the variation Srric = 0.2 GeV. The nPDF 
uncertainty is found by computing the asymmetric la 


^ These corrections were already implemented for similar predic¬ 
tions provided in other work Eoj. 
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CL from the eigenvectors basis, after scaling the eigen¬ 
vector deviations down by 1.645. The PDF uncertainty 
is found by computing the Icr CL from the replica set. 
A scale uncertainty, due to missing higher-order correc¬ 
tions, is evaluated by varying factorisation and renormal¬ 
isation scales independently by a factor of two around the 
nominal scale ^ with the constraint 1/2 < < 2. 

The nominal scale is set as /i = a/-I- p^, where px is 
the heavy-quark transverse momentum in the underly¬ 
ing Born conhguration. The sum of these contributions 
added in quadrature is taken to be the total uncertainty. 

Finally, it is worth noting that in very low-px region 
(< 1 GeV), extremely low scales are probed when per¬ 
forming scale variation. The baseline NNPDFS.OLHCb 
PDFs are only provided for scales above pp > ^ GeV, 
and the nPDFs pp) are only valid in the region 

{x > 10“'*, pp > 1.3 GeV), below which the (n)PDFs 
are frozen. For this reason, a lower pp cut of 1 GeV is 
enforced on observables which are p-r-integrated, which 
reduces the dependence of the observable on the poorly 
described low-Q^ region. An alternative approach would 
be to perform an extrapolation of (n)PDFs, which is now 
available for PDFs by default as of LHAPDF 6.1.5, and 
take the envelope of the two approaches as an uncer¬ 
tainty. 

CROSS SECTION PREDICTIONS AND 
KINEMATICS 

During the pPb run of the LHC in early 2013, LHCb 
recorded data samples corresponding to integrated lumi¬ 
nosities of 1.1 nb“* and 0.5 nb“* in forward and back¬ 
ward collisions respectively. The energy of the proton 
beam was measured to be Ep = 3988 GeV, while that 
of the Pb beam is Z/A ■ Ep = 1572 GeV per nucleon, 
resulting in a proton-nucleon GOM energy of = 

5008 GeV [37] (referred to as 5 TeV). Due to the asym¬ 
metry of the proton and nucleon beam energies in the lab 
frame, there is a rapidity shift between the GOM and lab 
frames of Ay = -1-0.465 for forward collisions and Ay = 
-0.465 for backward collisions. As the LHCb detector is 
only instrumented in the forward region 2.0 < rj < 4.5, 
this results in an effective coverage of 1.53 < y < 4.03 
and —2.47 < y < —4.97 for forward pPb collisions. For 
data taking during pPb collisions in Run-H, predictions 
for a proton beam energy of Ep = 6500 GeV, and corre¬ 
sponding Pb beam energy ot Z/A ■ Ep = 2563 GeV per 
nucleon, corresponding to a proton-nucleon GOM energy 
of ys/V = 8163 GeV have also been computed. These 
predictions are provided in the Appendix. 

Before providing double differential cross sections for 
D production, it is informative to first study how the 
forward acceptance and pp cuts effect the sampling of 
the nPDFs. In the lower region of Figure the nor¬ 
malised differential cross section with respect Log(xi) 


at 5 TeV is shown for low- and high-py regions. In both 
cases, the rapidity is limited to the COM overlap re¬ 
gion of forward and backward collisions (2.47 < ycOM < 
4.03). The blue lines indicate the sampling when a min¬ 
imum Pp cut of 4 GeV is required, which shifts both 
xi and X 2 to larger values as compared to the low-pp 
region — a consequence of producing a higher sys¬ 
tem. In the upper region of Figure |l] the gluon nPDF 
is shown as a function of x for values of Q = 3, 5 GeV 
which approximately correspond to the selected pp re¬ 
gions of the lower plot. In addition to that of EPS09 
(which is adopted in this work), it is also informative to 
consider the shape and magnitude of nPDFs provided by 
other groups. In the same Figure, the DSSZ [48] gluon 
nPDF modification at Q = 3 GeV is also shown for the 
region x > I0“*. The DSSZ central value is qualitatively 
similar to EPS09, in the sense that (anti-)shadowing is 
expected at low(moderate)-x values, with the main ex¬ 
ception that the magnitude of the nPDF corrections are 
substantially smaller. This behaviour is not supported by 
the forward J/i/ data [TSl[T3j, however, the main point 
is that D measurements at LHGb are sensitive to both 
the low-x {x < 10“^) shadowing regime, as well as the 
anti-shadowing regime for values of x ~ 0.05, and will 
therefore have large impact of future global analyses. It 
is also worth pointing out that the recent nGTEQ15 anal¬ 
ysis [33] prefers a stronger modification of the gluon PDF 
at low-x values than that of EPS09. The sensitivity of 
high-pT heavy quark measurements at LHGb to the gluon 
PDF at large-x has been discussed here [50] . 



FIG. 1: Upper region: NLO EPS09 gluon nuclear modifica¬ 
tion factor R^/pijiXg, Q) for Q = 3, 5 GeV. Lower region: Nor¬ 
malised differential cross section with respect to Log(xi). 
The is required to be within the COM overlap region of 
forward and backward collisions accessible to LHCb, and the 
effect of placing pp cuts on the is indicated. 

In the following, predictions are provided for dif¬ 
ferential cross sections in the Lab frame in forward colli¬ 
sions at 5 TeV. The fragmentation fraction /(c —>■ D°) = 
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0.565 130] is applied. Although only predictions are 
presented, it is reasonable to apply a ratio of fragmenta¬ 
tion fractions to obtain Ds, and predictions, as dif¬ 
ferences in the distribution of these different flavours due 
to fragmentation is negligible as compared to scale un¬ 
certainties. The results are presented in Figure]^ where 
the pt distributions for different rapidity bins are shown 
separately. Note that to visually distinguish these dif¬ 
ferent distributions, a rapidity dependent multiplication 
factor of 10“"* has been applied. 



FIG. 2: Differential D° cross section in forward pPb collisions 
at 5 TeV. Both the central value and total uncertainty for each 
rapidity bin is shown, and a multiplication factor of 10“"* has 
been applied to different rapidity bins as indicated on the plot. 

The uncertainties at \ow-pT are substantial, predomi¬ 
nantly due to large scale uncertainties arising at low , 
but also due to sizeable (n)PDF uncertainties. Conse¬ 
quently, direct constraints on the size of CNM effects in 
D production are unlikely to come from differential cross 
section measurements due to the dominance of scale un¬ 
certainties. Although only predictions for forward col¬ 
lisions are provided here at 5 TeV, double differential 
predictions including central values and total uncertain¬ 
ties for forward and backward collisions at 5 and 8 TeV 
are provided in the Appendix. 

RATIO PREDICTIONS 

A more suitable observable for assessing the size of 
CNM effects is Q — the ratio of cross sections in 
forward and backward collisions performed in the COM 
frame. This observable has the benefit that theoretical 
uncertainties partially cancel in this ratio, with the ex¬ 
ception of the nPDF uncertainties. The scale and charm 
mass uncertainties do not cancel exactly, since varying 
these parameters results in an altered sampling of the 
nPDFs. Predictions of Rfb will be presented differen¬ 


tially in Pt and y of the D hadron^. When theoreti¬ 
cal uncertainties are assessed for Rfb^ the variation of 
parameters is computed simultaneously for forward and 
backward differential cross sections. 



FIG. 3: Predictions of Rfb in pPb collisions with respect to ra¬ 
pidity for production at 5 TeV. Both the total uncertainty, 
and uncertainty excluding the dominant nPDF uncertainties 
are shown. The experimentally accessible region to LHCb is 
also highlighted. 

In Figure|^ differential Rfb predictions for hadrons 
with respect to rapidity at 5 TeV are provided. Both the 
total uncertainty, and uncertainty excluding the nPDF 
uncertainties are shown demonstrating that nPDFs un¬ 
certainties are entirely dominant. Across the accessible 
kinematic region of 2.47 < y < 4.03 (limited by the over¬ 
lap of forward and backward collisions), the central NLO 
EPS09 value ranges from 0.85 at low rapidity too 0.7 at 
high rapidity. At higher rapidity, the low-a; region of the 
nPDFs is sampled more frequently in forward collisions, 
and consequently this ratio is reduced due to shadow¬ 
ing. The LO prediction is systematically lower, primar¬ 
ily due to the fact that the NLO kinematics smear the 
nPDF sampling at low-a; to larger average values where 
the shadowing effect is less emphasised. 

In Figure similar predictions are provided differen¬ 
tially in Pt- The nPDF uncertainties are entirely dom¬ 
inant in this case too, however the central value is ob¬ 
served to be approximately flat with respect to pr- This 
behaviour can be understood by re-examining Figure [l] 
where the sampled region of the nPDFs for different px 
regions is shown. As the px cut of the 11° is increased, 
the sampled values of a;i, and X 2 increase. This par- 


® Given the high statistics expected in the available data set, mea¬ 
surements of Rfb performed double differentially in px and y 
may be feasible. Corresponding predictions are available on re¬ 
quest. 
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tially reduces the shadowing effect in forward collisions, 
but simultaneously increases the anti-shadowing effect 
in backward collisions. These two effects approximately 
compensate one another in the considered phase space 
region. The numerical values of these predictions are 
provided in the Appendix. 
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FIG. 4: The same as Figure but differentially in pr- This 
prediction is also for pPb collisions at 5 TeV. 


CONCLUSIONS 

Predictions for double differential cross sections and 
forward-backward ratios of D production in pPb colli¬ 
sions have been provided. These predictions demonstrate 
that D measurements at LHCb are sensitive to CNM ef¬ 
fects. In particular, it is shown that isolating specific pt 
regions within the LHCb acceptance can provide sensi¬ 
tivity to both shadowing and anti-shadowing. 

It has to be emphasised that the low-x region of 
nPDFs is currently not well constrained by data, and 
consequently different nPDF parameterisations can re¬ 
sult in conflicting predictions for the size CNM effects in 
the forward region. It is essential to perform measure¬ 
ments which extend the current constraints in the region 
X < 10“^. In addition to the proposed D measurements, 
low mass Drell-Yan measurements should also be per¬ 
formed [Sn 152] . 
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APPENDIX 


The double differential cross sections for forward and 
backward production in the Laboratory frame at 5 
and 8 TeV are provided in the following Tables. The NLO 
predictions for Rft provided in Figure and Figure]^ 
are also provided, and in addition those at 8 TeV. As 
previously mentioned, double differential predictions of 
Rfb are available upon request. 


A >Pl^) A / L //-I \T\ 

^ - ...oAd (mb/GeV) 


P ? (GeV) 



2.0 

- 2.5 

2.5 

- 3.0 

3.0 

- 3.5 

3.5 

- 4.0 

4.0 

- 4.5 

0.0 

- 1.0 

6.43 

+ 15.61 
- 5.75 

6.11 

+ 15.04 
- 5.31 

5.86 

+ 14.55 
- 5.02 

5.37 

+ 13.51 
- 4.51 

4.67 

+ 11.87 
- 3.89 

1.0 

- 2.0 

9.33 

+ 19.76 
- 7.99 

8.95 

+ 19.08 
- 7.55 

8.27 

+ 17.93 
- 6.73 

7.43 

+ 16.29 
- 5.91 

6.48 

+ 14.58 
- 5.11 

2.0 

- 3.0 

6.2 

+ 10.85 
- 4.95 

5.81 

+ 10.29 
- 4.52 

5.24 

+ 9.41 
- 3.98 

4.54 

+ 8.11 
- 3.43 

3.74 

+ 6.92 

- 2.76 

3.0 

- 4.0 

3.39 

+ 4.98 

- 2.46 

3.11 

+ 4.65 

- 2.22 

2.74 

+ 4.12 

- 1.92 

2.28 

+ 3.49 
- 1.55 

1.76 

+ 2.77 

- 1.16 

4.0 

- 5.0 

1.79 

+ 2.33 
- 1.19 

1.61 

+ 2.12 
- 1.05 

1.37 

+ 1.81 
- 0.87 

1.1 

+ 1.51 
- 0.68 

0.81 

+ 1.07 
- 0.49 

5.0 

- 6.0 

0.95 

+ 1.08 
- 0.58 

0.82 

+ 0.94 
- 0.48 

0.67 

+ 0.79 
- 0.38 

0.53 

+ 0.64 
- 0.3 

0.39 

+ 0.47 

- 0.22 

6.0 

- 7.0 

0.51 

+ 0.54 

- 0.28 

0.46 

+ 0.49 
- 0.25 

0.35 

+ 0.39 
- 0.18 

0.27 

+ 0.3 
- 0.14 

0.19 

+ 0.21 

- 0.1 

7.0 

- 8.0 

0.3 

+ 0.29 

- 0.15 

0.27 

+ 0.25 
- 0.13 

0.2 

+ 0.19 
- 0.1 

0.15 

+ 0.15 
- 0.07 

0.09 

+ 0.09 

- 0.04 

8.0 

- 9.0 

0.17 

+ 0.16 

- 0.08 

0.15 

+ 0.14 
- 0.07 

0.12 

+ 0.11 
- 0.05 

0.09 

+ 0.08 
- 0.04 

0.05 

+ 0.05 

- 0.02 

9.0 

- 10.0 

0.11 

+ 0.09 

- 0.05 

0.1 

+ 0.08 
- 0.04 

0.07 

+ 0.07 
- 0.03 

0.05 

+ 0.05 
- 0.02 

0.02 

+ 0.02 

- 0.01 


TABLE I: Predictions for the differential cross-sections for 
meson production in forward pPb collisions at 5 TeV. 


(r-b/GeV) 


P? (GeV) 





y 

D 






2.0 

- 2.5 

2.5 

- 3.0 

3.0 

- 3.5 

3.5 

- 4.0 

4.0 

- 4.5 

0.0 

- 1.0 

6.86 

+ 16.56 
- 6.16 

6.83 

+ 16.69 
- 6.07 

6.74 

+ 16.73 
- 5.87 

6.59 

+ 16.71 
- 5.6 

6.25 

+ 16.11 
- 5.31 

1.0 

- 2.0 

10.02 

+ 21.0 
- 8.66 

9.91 

+ 21.08 
- 8.36 

9.82 

+ 21.17 

- 8.13 

9.33 

+ 20.49 
- 7.58 

8.63 

+ 19.44 
- 6,87 

2.0 

- 3.0 

6.7 

+ 11.54 
- 5.41 

6.58 

+ 11.55 
- 5.17 

6.3 

+ 11.21 
- 4.83 

5.83 

+ 10.6 
- 4.37 

5.09 

+ 9.48 

- 3.76 

3.0 

- 4.0 

3.77 

+ 5.58 
- 2.77 

3.55 

+ 5.29 
- 2.54 

3.32 

+ 4.96 

- 2.32 

2.97 

+ 4.59 
- 2.01 

2.43 

+ 3.9 

- 1.64 

4.0 

- 5.0 

1.96 

+ 2.48 
- 1.3 

1.86 

+ 2.4 
- 1.21 

1.68 

+ 2.19 

- 1.05 

1.42 

+ 1.9 
- 0.87 

1.12 

+ 1.53 
- 0.67 

5.0 

- 6.0 

1.06 

+ 1.21 
- 0.64 

0.98 

+ 1.13 
- 0.58 

0.84 

+ 0.97 

- 0.48 

0.68 

+ 0.81 
- 0.37 

0.49 

+ 0.61 

- 0.26 

6.0 

- 7.0 

0.57 

+ 0.6 
- 0.32 

0.53 

+ 0.57 
- 0.28 

0.44 

+ 0.48 

- 0.23 

0.36 

+ 0.39 
- 0.18 

0.24 

+ 0.27 

- 0.12 

7.0 

- 8.0 

0.34 

+ 0.31 
- 0.17 

0.31 

+ 0.3 
- 0.15 

0.25 

+ 0.25 

- 0.12 

0.19 

+ 0.19 

- 0.08 

0.11 

+ 0.12 

- 0,05 

8.0 

- 9.0 

0.21 

+ 0.19 
- 0.09 

0.18 

+ 0.17 
- 0.08 

0.15 

+ 0.15 

- 0.07 

0.1 

+ 0.1 

- 0.04 

0.06 

+ 0.05 

- 0.02 

9.0 

- 10.0 

0.13 

+ 0.11 
- 0.06 

0.11 

+ 0,1 
- 0.05 

0.09 

+ 0.08 

- 0.04 

0.06 

+ 0.06 

- 0.02 

0.03 

+ 0.03 

- 0.01 


TABLE H: Predictions for the differential cross-sections for 
meson production in backward pPb collisions at 5 TeV. 



dRfbldx 

P? (GeV) 

1.5 

3.0 

5.0 

7.0 

9.0 

NLO EPS09 

n Cl + 0.06 

0-81 o 06 

0 79 

^ - 0.06 

o 

oo 

1 + 
O o 

b b 

-81 -1 

0 79 

‘ ^ - 0,07 

0 82 

- 0.08 

y 

2.25 

2.75 

3.25 

3.75 

4.25 

NLO EPS09 

n so 

0.89 _ o . o 3 

n oA + 0.06 

U - o 4 _ 0,05 

0 80 

U.OU _o 06 

0 75 

^ - 0.07 

0 71 + 0'09 
^ - 0.09 


TABLE HI: Differential predictions for Rfb in pPb collisions 
at 5 TeV. 
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TABLE IV: Predictions for the differential cross-sections for 
_D° meson production in forward pPb collisions at 8 TeV. 


(^fc/GeV) 



TABLE V: Predictions for the differential cross-sections for 
_D° meson production in backward pPb collisions at 8 TeV. 



dRfhjdx 

Pt (GeV) 

1.5 

3.0 

5.0 

7.0 

9.0 

NLO EPS09 

n +0.06 

-0.05 

n 8R +0.06 
0.33 _0,05 

o 

bo 

to 

1 + 
o o 
b b 

0 79 

‘ ^ -0.06 

n 78 +0.07 
O-'o -0.07 

y 

2.25 

2.75 

3.25 

3.75 

4.25 

NLO EPS09 

0 92 

-0.02 

o 

bo 

GO 

1 + 
O o 

b b 

CO 

o o 
b b 
+ 1 

+ 

oo 

6 

0 79 

‘ ^ -0.06 

0 74 +0-07 
-0.07 


TABLE VI: Differential predictions for Rfb in pPb collisions 
at 8 TeV. 



